New interpretation of the origin of 2DEG states at the surface of layered topological 
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On the basis of relativistic ab-initio calculations we show that the driving mechanism of simulta- 
neous emergence of parabolic and M-shaped 2D electron gas (2DEG) bands at the surface of layered 
topological insulators as well as Rashba-splitting of the former states is an expansion of van der 
Waals (vdW) spacings caused by intercalation of metal atoms or residual gases. The expansion 
of vdW spacings and emergence of the 2DEG states localized in the (sub)surface region are also 
accompanied by a relocation of the topological surface state to the lower quintuple layers, that can 
explain the absence of interband scattering found experimentally. 

PACS numbers: 73.20.-r, 79.60.-i 
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The recently discovered three-dimensional topological 
insulators belong to a class of insulators in which the 
bulk gap is inverted due to the strong spin-orbit interac- 
tion [l|-|j] . A direct consequence of such bulk band struc- 
ture arises at the surface: the spin-polarized topologically 
protected massless metallic states, forming a Dirac cone. 
These surface states (SS) exhibit many interesting prop- 
erties resulting from the fact that the spin of electron is 
locked perpendicular to its momentum, thus forming a 
SS spin structure that protects electrons from backscat- 
tcring under the influence of weak perturbations. This 
makes topological insulators potentially promising mate- 
rials for creation of new quantum devices. 

By now few families of topological insulators (TI) have 
been discovered [5H15| , of which the binary layered com- 
pounds Bi2Te3, BJ2Se3, and Sb2Te3 are the most studied 
both experimentally [l2-[l5| and theoretically [1, [H-IIoj] . 
The latter systems have tetradymitc-like layered struc- 
ture with ionic-covalent bonded quintuple layer (QL) 
slabs, which are linked by weak van der Waals forces. 
Such a layered structure predetermines the formation of 
the surface by cleavage on the van der Waals (vdW) 
spacing that doesn't result in the formation of dangling 
bonds, and thus only the Dirac states reside in the bulk 
energy gap. 

However, recently in several studies it has been demon- 
strated by using angle-resolved photoemission spec- 
troscopy (ARPES) that beside the Dirac cone the 2DEG 
states arise at the surface of Bi 2 Se3 after a few hours 
of exposition in vacuum [2l| , upon deposition of various 
magnetic 23J, [24 1 and non-magnetic atoms 241426 1. and 
molecules as well 25|, |27j . These states form a parabolic 
band (PB) in the energy gap just below the conduc- 
tion band and M-shapcd band in the local gap of bulk- 



projected valence band. The former bands show an ap- 
preciable Rashba spin splitting [22J42J, |26j, |27 1 . Moreover, 
for several adsorbates 2J,l26j at the saturation deposition 
time a second and even third pair of spin-split parabolic 
states emerges below the conduction band minimum. In 
the most of these papers the emergence of PB states was 
ascribed to a confinement of the conduction band states 
in a quantum well formed by band bending potential pro- 
duced by adsorption of metallic atoms or residual gases, 
although in Ref. 27 it was pointed out that the poten- 
tial gradient from band bending cannot alone be respon- 
sible for their Rashba-splitting. Indeed, a model calcu- 
lation exploiting band bending approach based on the 
coupled solution of the Poisson and Schrodinger equa- 
tions yielded the Rashba spin-orbit coupling parameter 
an which is significantly smaller than that extracted from 
the ARPES data [22|. Moreover this approach doesn't 
reproduce the M-shaped band [21] . The most surprising 
experimental finding is the absence of interband scatter- 
ing of the Dirac state electrons in the presence of the 
PB spin-split states after the deposition of any kind of 
metal atoms |2J]. A study of a naturally aged surface 
of Bi2Se3 also reveals that scattering rate in the topo- 
logically protected state is unaffected by the potential 
created by adsorbed atoms or molecules 28 1 . 



In this Letter we propose an alternative explanation of 
the emergence of the 2DEG states based on a well known 
fact: the intcrlayer gaps (vdW spacings) in the layered 
compounds can serve as natural containers for impurities 
in synthesis processes and for intercalated atoms. Vari- 
ous atoms have been intercalated in different layered ma- 
terials, inducing expansion of vdW spacings 29143 If. One 



would expect that due to a weak binding between QLs 
even a relatively small concentration of contaminants in 
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FIG. 1: (Color online) (a) Surface band structure of Bi 2 Se3 
in the vicinity of f with the vdW gap expanded from 10% 
to 50% (insert: variation of the lower part of the Dirac cone 
with vdW gap expansion); (b) a magnified view of Rashba- 
split parabolic bands. 



the vdW gap can produce its sizeable expansion. The 
magnitude of this expansion depends on the impurity 
size and on formation of impurity atom clusters within 
the vdW gap. Mechanisms of diffision of adsorbed atoms 
into the vdW gap are not known in detail yet, however, 
at least two kinds of the diffusion process can be antici- 
pated. First, the Fermi level of naturally grown crystals 
of the tetradymite-like TI is usually found to be located in 
the bulk conduction band due to vacancy defects [l2|, Il5| 
and, thus, apparently there is a sufficiently high vacancy 
concentration to ensure the vacancy mediated penetra- 
tion of impurities through a QL into the vdW gap. The 
second process can be associated with the sliding of im- 
purity atoms into the vdW gap from step edges. 

Here we show on the basis of ab-initio calculations us- 
ing the VASP code [32, 



that the driving mechanism 
of the simultaneous emergence of the parabolic and M- 
shaped bands in Bi2Sc3, Bi2Te3, and Sb2Te3, is a widen- 
ing of the outermost vdW spacing. This scenario also 
explains the Rashba-type splitting for the PB. We find 
that besides the emergence of PB states, which are lo- 
calized in the detached QL the expansion of the vdW 
spacing results in a relocation of the Dirac state to the 
lower QL that makes topological and PB states sepa- 
rated in space. This fact explains the observable absence 

" r2? 



of interband scattering [24 



We also show that the 
expansion of various vdW gaps produces multiple 2DEG 
states. 

We start with consideration of the Bi2Se3 surface. To 
simulate the effect of the adsorbate deposition time as 
well as the impurity atom size we perform a calculation of 



TABLE I: The energy with respect to Ef, the fitted effec- 
tive mass, and the Rashba coupling parameter ota for the 
parabolic SS as a function of the vdW gap expansion. 



expansion 
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Bi 2 Se 3 



10 
20 
30 
40 
50 



0.30 
0.27 
0.24 
0.22 
0.19 



0.38 
0.27 
0.20 
0.18 
0.17 



0.46 
0.28 
0.20 
0.18 
0.17 



0.16 
0.25 
0.33 
0.29 
0.24 



0.14 
0.24 
0.32 
0.29 
0.24 
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10 


0.15 


0.31 


0.37 


0.29 


0.28 




20 


0.04 


0.17 


0.18 


0.40 


0.40 




30 


0.01 


0.14 


0.13 


0.74 


0.71 




40 


0.00 


0.12 


0.12 


0.59 


0.57 




50 


-0.02 


0.12 


0.13 


0.45 


0.45 



Sb 2 Te 3 



10 
20 
30 

40 
50 



0.09 
0.08 
0.05 
0.04 
0.02 



0.24 
0.17 
0.12 
0.12 
0.12 



0.28 
0.18 
0.13 
0.12 
0.12 



0.30 
0.35 
0.26 
0.15 
0.10 



0.30 
0.36 
0.26 
0.15 
0.10 



the Bi2Se3 surface with expansion of the outermost vdW 
spacing by 10% - 50%. As one can see (Fig. Q]), the de- 
tachment of the outermost QL leads to the simultaneous 
emergence of both the Rashba-split band below the bot- 
tom of the bulk conduction band and the M-shaped band 
in the valence bulk-projected gap. The energy of the 
M-shaped and parabolic spin-split bands as well as the 
magnitude of the Rashba splitting parameter strongly de- 
pend on the vdW spacing expansion. At 10% expansion 
the first pair of spin-split bands as well as the M-shaped 
state emerges just below the bulk projected states. Upon 
increasing the vdW expansion these bands shift gradu- 
ally down. This behavior reflects the rise of the 2DEG 
bands as a function of the deposition time found in re- 
cent experimental studies 2J,|27|. Up to 30% expansion 
this also leads to the decrease of the effective mass m* of 
PBs and to the increase of the Rashba coupling parame- 
ter an (see Table H| . At the same time m* demonstrates 
an apparent L-K/L-M anisotropy at 10% expansion. At 
higher expansions this anisotropy is significantly reduced. 
The small (« 0.01) m* anisotropy was also obtained in 
Ref. [Hi (SI). Different fitting procedures for the ARPES 
measured band in Ref. |22| (SI) gave the Rashba splitting- 
parameter varying from 0.36 to 1.35 eV-A in aged Bi 2 Sc 3 
surface. In the case of potassium deposition on Bi2Se3 a 
maximal a,R — 0.79 ± 0.03 eV-A was observed after 2.5-3 
minute K evaporation [26j . Our model gives a lower limit 
for the former case at moderate expansions and underes- 
timates the latter maximal an while the band bending 
model 22| gives 0.1 eV-A for the PB splitting. 

The nature of spin-split PBs can be deduced from anal- 
ysis of their spacial localization. At 10% expansion the 
emerged parabolic state, which splits off from the con- 
duction bulk band continuum, keeps yet a bulk-like delo- 
calized character (Fig. [^a)). This explains the obtained 
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FIG. 2: (Color online) (a), Spatial localization of the topolog- 
ical (black) and parabolic (red) surface states as a function of 
the vdW gap expansion; green line represents localization of 
the M-shaped band for 50% vdW gap expansion; (b), charge 
density distribution of the PB state for 30% expansion; (c) 
and (d), |^(z)| 2 of the lowest unoccupied band for a free- 
standing 1-QL and of the upper occupied band for a free- 
standing 2-QL with 50% vdW gap expanded, respectively. 



f-K/f-M anisotropy of the parabolic band at small ex- 
pansions: it just reflects a similar anisotropy of the bulk 
states. Upon increasing the vdW spacing this state ac- 
quires more localized character and at 30-50% it almost 
completely lies in the detached QL. The PB state, char- 
acterized by four peaks localized in the outermost QL, 
Fig. [Ha), reflects the symmetry of the energy gap edges 
at T: it is composed of the Bi p z orbitals (lower edge 
of the gap) with the sizeable Se p z contribution (upper 
edge of the gap) to the vdW-side peak (Fig. EJb)). With 
the further increase of the outermost vdW spacing the 
localization of the parabolic state approaches the distri- 
bution ^(z)! 2 of the lowest unoccupied band in a free- 
standing QL (Fig.[2{c)). The latter band is not spin-orbit 
split since a free-standing QL possess inversion symme- 
try. This effect accounts for the reduction of an in the 
PB states with the increase of the vdW spacing beyond 
30% (Tabic HI . Note here that a large concentration of 
intercalated atoms in the vdW gap can affect the vdW- 
side Se orbitals and hence lead to change of some features 
of the spin-split band, for instance increase the Rashba 
splitting. 

In the case of Bi2Te3 and Sb2Te3 we find the same 
trend in the development of the Rashba-split PB (Ta- 
ble H]). However, the heavier Te provides two times big- 
ger an in Bi2Te3 with respect to that in Bi2Se3 while 
the lighter Sb, in turn, gives a smaller Rashba coupling 
parameter in Sb2Tc3. 
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FIG. 3: (Color online) Surface band structure of Bi2Se3 with 
the first and second vdW gaps expanded by 40% and 10%, 
respectively (a); by 40% and 20%, respectively (b); spatial 
localization of the topological state (black line) as well as for 
the lower (red) and upper (blue) Rashba-split SSs for the (b) 
case (c). 



An intriguing finding is that the development of PB 
states is accompanied by shifting the topological state 
deep into the crystal, so that at 30-50% expansion the 
Dirac state is mostly located in the second QL, beneath 
the detached QL. In this case wave functions of the topo- 
logical and parabolic states only slightly overlap. The 
latter explains the experimentally observed absence of 
interband scattering 2J, |28[ . It should be stressed that 
the dispersion of the upper part of the Dirac cone as well 
as the position of the Dirac point remain unchanged un- 
der this relocation while the dispersion of the lower part 
of the Dirac cone is slightly modified (Fig. [TJa), insert) 
in agreement with ARPES data [25 1 . 

The M-shaped band emerges in the local gap of the 
bulk-projected valence band splitting off from the up- 
per edge of the gap. The edge is formed by the bulk 
Se p z states only which arc not spin-orbit split at T and 
around. The M-shapcd SS being split off from the Se p z 
bands does not acquire other symmetry and thus does not 
show spin-orbit splitting at small momenta. Only being 
close to the bulk continuum states (at large momenta) it 
shows weak spin-orbit splitting. With the increase of the 
vdW gap it becomes more localized. Charge density of 
this state as shown for 50% expansion case in Fig. [Ha) is 
mostly situated in the expanded vdW gap and around. 
Comparing this ^(z)! 2 with that of the the upper oc- 
cupied state in free-standing 2-QL with 50% vdW gap 
expanded, Fig. WiA), one can note that these two states 
are of the same origin: the vdW spacing expansion. 

Now we consider the effect of simultaneous expansion 
of the first and second vdW spacings. The expansion of 
the latter can be caused by further diffusion in the sample 
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FIG. 4: (Color online) Spin structure of the Dirac state and 
Rashba-split states of Bi2Se3, as given by spin projections S x , 
S y , and S z at energy of 280 meV for the detached 1QL (30%) 
case (a) and for the 2QL (40% + 20%) case (b). 



of a smaller amount of atoms deposited on the surface. 
In Fig. [3Ja) we show the band spectrum of Bi2Se3 with 
40% and 10% expanded 1-st and 2-nd vdW spacings, 
respectively. This film geometry gives rise to the second 
parabolic state below the conduction band. This state 
shows very small spin splitting. The emergence of the 
second PB leads to a decrease of energy of the first band 
by 20 meV as well as a reduction of its atR to 0.27 cV-A 
compared to the case of the single vdW spacing expansion 
by 40%. Further increase of the 2-nd vdW gap leads to 
a subsequent downshift of both PBs (Fig. Etb)). While 
an for the upper band in this case is 0.16 (0.15) eVA in 
f-K (f-M) directions, respectively, the spin splitting of 
the lower state is reduced to 0.25 eV-A. Localizations of 
the lower and upper bands in general resemble those for 
a single vdW expansion by 40% and 20%, respectively 
except their mutual overlap (Fig. EJc)). The overlap of 
the topological SS with that localized in the 2-nd QL is 
similar to that in the single 20% case while the Dirac 
state is absent in the outermost QL as in the case of 
the single 40% detachment. Further increase of the 2- 
nd vdW spacing leads to bigger downshift of both PBs, 
to increase of their overlap and complete relocation of 
the Dirac state to the third QL. The expansion of the 
second vdW spacing also results in the development of 
the second M-shaped state in the valence band gap (not 
shown). Subsequently expansion of the third vdW gap, 
that simulate a deeper diffusion of deposited atoms, leads 
to simultaneous emergence of the third pair of parabolic 
Rashba-split and M-shaped bands (not shown). 

Note that the sizeable overlap between the topologi- 
cal SS and 2DEG state exists barely at the beginning of 
the formation of each new PB and hence the interband 
scattering is possible in a narrow energy interval just be- 



low the conduction band only. In fact this scattering 
will be also limited by spin conservation rule due to spin 
polarized nature of both topological and 2DEG Rashba- 
split states. The Dirac state has clockwise spin helicity 
while the parabolic states demonstrate typical Rashba- 
type counter-clockwise and clockwise helicities for outer 
and inner branches, respectively, with a very small S z 
spin component for both types of SS's (Fig. HJa)). This 
spin structure allows "ki to k2" scattering between the 
topological state and the inner branch of the PB and "ki 
to — k2" transitions between the Dirac cone state and 
the outer branch of the PB. However, the efficiency of 
these transitions should be significantly reduced because 
of small overlap between the topological SS and PB state 
for well developed PBs. In the case of multiple detached 
vdW gaps new channels for scattering of Dirac electrons 
(Fig. Stb)) arise, nevertheless, only PBs that have size- 
able overlap with the topological SS can contribute to 
considerable interband scattering. 

To summarize, on the base of the proposed scenario 
of impurity intercalation in the vdW gaps of the layered 
TIs and ab-initio calculation results we have shown that 
the driving mechanism of experimentally observed simul- 
taneous emergence of parabolic and M-shaped bands as 
well as Rashba-splitting of former bands is expansion of 
vdW spacings. These bands are 2DEG states localized in 
the narrow (sub)surface region. We have found that be- 
sides the development of 2DEG SS's the expansion of the 
vdW gap also provides relocation of the topological state 
to the lower QL. This fact explains the observed absence 
of interband scattering between the topological SS and 
Rashba-split PB states. We also show that expansion of 
the low-lying vdW gaps produces multiple 2DEG states. 
The proposed mechanism is expected to be also valid for 
other layered TIs. 
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